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Abstract 

Effects of geometric constraints on a steady flow potential are described 
by an elliptic-hyperbolic generalization of the harmonic map equations. 
Sufficient conditions are given for global triviality. MSC2000: 58E20, 
58E99, 75N10. 



1 Introduction 

In |l(Jj , local properties of maps which are critical points of a nonlinear Hodge 
energy were investigated (see also [8]). The target of the map has a physical 
interpretation as a geometric constraint on the potential of a steady flow. In 
this note we illustrate these maps by deriving elementary but explicit physical 
examples, clarify the relation of the objects studied in ^H] to other classes of 
maps which have been studied recently, and provide global conditions under 
which critical points reduce to the trivial map. 

1.1 A column of tap water as a mapping 

A simple motivating example for placing geometric constraints on a flow poten- 
tial is provided by the steady flow of water from a faucet. If Vi is the velocity 
of the flow through a thin horizontal section of area Ai and if Vf and Af are 
defined analogously, then the conservation of mass implies that 



ViAi=VfA f . (1) 

But the particles accelerate under gravity, so if the cross section Af is taken 
nearer to the drain than the cross section Ai, we conclude that Vf > Vi. Equation 
l|T]l then implies that Af < Ai, which explains why the column of water is seen to 
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taper. In this conventional approach the flow geometry is derived by imposing 
a physical law (acceleration under gravity) on a conservation law. 

An alternative approach would be to observe that the column of water tapers, 
so that Af < Ai. Conservation of mass implies eq. (I), so we conclude that 
Vf > Vi, that is, the particles accelerate under gravity. In this alternative 
approach the physical law is derived by imposing an observed flow geometry on 
a conservation law. 

It is convenient to think of the geometry as a constraint applied to a flow 
potential u. In this example u maps a right circular cylinder into a tapered 
cylinder. Recall that every smooth curve has a dual representation as the enve- 
lope of its family of tangent lines. In a steady flow the velocity vectors appear 
as tangent lines to the potential surfaces, and so form an envelope of the cross 
sections represented in eq. QJ. We see these cross sections taper by tracing 
the velocity vectors of the water droplets. Thus, as is often the case in fluid 
dynamics, the mathematical abstraction of a potential surface attains a visible 
representation in physical space. 

While this simple example involves incompressible flow, the same alterna- 
tives exist in the more complicated case of compressible flow such as the flow of 
exhaust from a jet engine. 

1.2 Shallow hydrodynamic flow 

Now we consider a slightly more complicated case, that of steady, inviscid, hy- 
drodynamic flow in a shallow channel. Write the flow velocity v in components 

t>2, t>3) , where v± is the horizontal component in the x-direction, vi is the 
horizontal component in the y-direction, and V3 is the component in the (verti- 
cal) z-direction. Impose initial conditions under which 1)3 is zero at time t = 0. 
Because we are assuming shallow depth, it is reasonable to suppose that the 
component of acceleration of water particles in the z-direction has negligible ef- 
fect on pressure. The result of applying this hydrostatic law is that 1*3 remains 
zero for all subsequent times and the horizontal velocity components V\ and 1*2 
are independent of the z-coordinate. Because the flow is steady, the velocity 
components are also independent of t. 

Generalizing eq. to express the vanishing of an appropriate surface inte- 
gral and applying the Divergence Theorem, we write the law of mass conserva- 
tion in the form of a continuity equation (see, e.g., |HJ, Sec. 1.1.1) 

d d 

— [h (x, y) vi (x, y)\ + — [h (x, y) v 2 (x, y)] = 0, (2) 

where h(x, y) represents the depth of the channel at the point (x, y). Bernoulli's 
formula expresses h as a function of Q = \v\ 2 , that is, h (Q) — (C — Q) /2g, 
where C is a constant and g is the magnitude of gravitational acceleration. 
Substituting this relation into J5J and using the chain rule, we find that prior 
to the imposition of any geometric constraint the flow will satisfy 
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C-Q 
2 



v\ V lx - V X V 2 (Vly + V 2x ) + 



C-Q 
2 




= 0. 



(3) 



(In eq. JHJ numerical subscripts denote vector components, whereas variable 
subscripts denote partial differentiation in the direction of the variable.) When 
we applied Bernoulli's formula we tacitly assumed the flow to be irrotational. 
Thus its velocity vector has vanishing curl, which allows us to equate mixed 
partial derivatives and assume the local existence of a potential function u(x, y) 
such that Vu = v. 



Writing c 2 = gh = (C-Q) /2, we obtain (c.f. (10.12.5) of Q7|) a second- 
order quasilinear elliptic-hyperbolic equation for the potential function: 



The type of eq. Q depends on whether or not the flow speed \[Q exceeds the 
propagation speed c. For subcritical flow speeds in which the Froude number 
F = \fQ /c is exceeded by 1, the continuity equation is of elliptic type and the 
flow is tranquil. For supercritical flow speeds in which F exceeds 1, eq. Q is of 
hyperbolic type, which characterizes shooting flow. 

We can prescribe flow geometry for this problem by means analogous to 
the simpler example of water flowing from a faucet. In this case the potential 
function can be considered as a map from the flow domain to a target manifold. 
In order for the map to have non-trivial geometry, we must assume that the 
potential function it is multi-valued. Its components can be imagined as local 
coordinates on the target manifold. However, unlike the usual representations 
of the flow potential in the complex plane (e.g., Sec. 1.6 of 6 ), our mappings 
are defined over the real field. A geometric variational problem for this model 
will be formulated in the next section. Note that multi-valued potentials arise 
naturally on flow domains having non-trivial topology. 

Conjectures on the observable effects of geometric constraints on the flow of 
shallow water go back at least 400 years. Geometric arguments, to one degree 
or another, have been applied to explain tidal bores on the English rivers Severn 
and Trent, the French river Seine near Caudebec-en-Caux, and the Chinese river 
Tsien-Tang, as well non-tidal anomalies such as those involving the Agulhas 
Current. See |12j for a review. As our last physical example, we recall a model 
for those effects which has particularly simple geometry; the model is given in 
greater detail in |18j . 

Consider a steady current flowing in the positive- a; direction. Suppose that 
an incline of magnitude 8 occurring between the points xq and x\ of an otherwise 
horizontal channel floor produces a surface elevation of height e at x = X\. 
Suppose that xo < X\, that the velocity of the flow to the left of xq is v\ 
and that the velocity of the flow to the right of x\ is v\. An inviscid flow is a 
conservative system. Equating the kinetic and potential energies of the flow for 
an arbitrary surface particle of mass m, we have 





U 2 y ] Uyy = 0. 



(4) 



mvf mvi 



(5) 
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Write the continuity equation for this one-dimensional system in the form 



{H + S)v 1 = (H + e)v 1 . (6) 
Equations and JBJ can be combined into the single expression 

(H 2 + 2He + e 2 ) 2ge = v\ (2H + e + 6) (e - 8) . 

Approximating this expression to first order in e and <5, we obtain 

- 6 - 6 

1 - gH/vl ~ l-(c/ Vl f 

Again we find that the character of the flow is determined by whether the Froude 
number F exceeds, equals, or is exceeded by the number I , where in this case 
y/Q = \vi | . The blow-up singularity that develops as F tends to f is avoided by 
hypothesis, as e is assumed to be small. In the case of tranquil flow, a positive 
elevation |e| occurs when S is exceeded by zero and a negative elevation — |e| 
occurs when S exceeds zero. The opposite relations hold for shooting flow. Of 
course the effects of turbulence are ignored. 



2 A geometric variational problem 

Equation (@J can be derived by a variational principle from an energy functional 
having the form 

E p (u;Z,R n )= [ [ Q p(s)dsdE, (7) 

JT, JO 

where E is a surface and p(Q) = c 2 . In generalizing this problem we consider 
an energy functional of the form 

r rQ(du) 

E p (u; M, N) = I I p(s) dsdM, (8) 

J M JO 

where M is a Riemannian manifold of dimension n; N is a Riemannian manifold 
of dimension m; u : M — > N is a bounded map; 

Q {du) = (du, du) | T , Mm -i TN ; (9) 
p : K + U {0} — > K + is a C 1,a function of Q satisfying the differential inequality 

0< M^l <m (10) 

p(Q) 

for Q 6 [0, Q C rit)', Qcrit is the square of the sonic flow speed. Inequality (|I0|l is 
a condition for tranquil channel flow. 

We are interested in maps u which are stationary with respect to E p (u; M, N) 
in the sense that 



4 



5E = 



-E(u t ;M,N)\ t=0 = 0, 



(11) 



where u t : M x (— e < t < e) — > iV, u = u, is a smooth, compactly supported 
one-parameter family of variations (to be further specified below). 

The variational equations of E p (it; M, N) are satisfied by maps which are 
extremal within a competing homotopy class of finite-energy maps from M to 
N. However, there are solutions of which are not extremal with respect to 
any class of maps. 

Condition IjlOjl is a condition for ellipticity of the variational equation for E, 



where V co „ denotes the covariant derivative in the bundle T*M(£>u 1 TN. Note 
that this equation introduces geometry into both the domain and range of eq. 
0- 

As an alternative to the hydrodynamic interpretation, the manifold N may 
be chosen to represent a geometric constraint placed on the flow potential of a 
steady, irrotational, polytropic, perfect, compressible fluid, which is adiabatic 
and isentropic and for which the closed 1-form rfueT (T*M) is dual to the flow 
velocity. In that case we choose 



where 7 > 1 is the adiabatic constant of the medium [H] . Those choices trans- 
form l|lUfl into a condition for subsonic compressible flow of mass density p. The 
sonic transition as Q tends to Q cr it is a gas-dynamic analogy for the change in 
the aspect of tap water from clear to white at a sufficiently high velocity or for 
the transition from tranquil to shooting channel flow in hydraulics. We recover 
harmonic maps as the incompressible limit p{Q) = 1. 

Note that two distinct terms are both referred to as density in the mathe- 
matical/fluid dynamics literature. The physical, or mass density of the flow is 
p(Q), but the density of the variational integral E p is given by 



In particular, the physical density given by l|13|l is a decreasing function of 
Q, whereas the corresponding variational density is an increasing function of 
Q provided Q < 2/(7 — 1). The variational density e(u) corresponding to the 
physical density (fC^ vanishes (or cavitates) at the flow speed Q — 0, but the 
physical density itself does not. The physical density Ijl3(l cavitates at the flow 
speed Q = 2/(7 — 1), but the variational density does not. 

Although in realistic physical contexts the ratio 7 of specific heats must be 
taken to exceed 1, many of the analytic properties of formula (|13|l extend to 



traceV cov (p(Q)du) = 0, 



(12) 




(13) 
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the limiting case in which 7 tends to unity. In this limiting case the variational 
density is given by the function 

-e{u) = 1 - exp [-Q/2] . 

(We can easily see this by writing y = \ogp and applying L'Hopital's rule to 
QRijl.) Regularity arguments for weak subsonic solutions of presented for the 
unconstrained case in, e.g., can be extended to the limiting value of 7 by the 
Arzela-Ascoli Theorem. One thus obtains the existence of continuous subsonic 
solutions in that limit by semicontinuity, using the convexity of the limiting 
energy on the subsonic range (c.f. Proposition 1 of [HI)- Some comments on the 
geometry of the limiting case are given in Sec. 2.1 of jH]. 

Reference includes a review the isometric embedding of the target man- 
ifold for critical points of JSJ, © into a higher-dimensional Euclidean space K fe 
and the use of nearest-point projection to obtain a form of the geometric con- 
straint which is convenient for variational analysis in a Sobolev space. This is a 
familiar trick in the theory of harmonic maps j!4| . We obtain variations of the 
form ut — 7Tjv {u + tip) , where ifi is a smooth map from M into M. k and 7rjy is 
the nearest-point projection, assigning to every y in a Euclidean neighborhood 
of N the point on TV that minimizes the distance to y. The embedding of the 
flow geometry in Euclidean space does not necessarily embed the corresponding 
physics in an ambient Euclidean space. To illustrate the distinction, compare 
water poured over a small sphere on the surface of Earth with water flowing on 
the surface of a spherical planet. In the former case we would take the gravi- 
tational acceleration vector to point in the vertical direction of the Euclidean 
coordinate system in which the sphere sits; in the latter case we would take the 
gravitational acceleration to point in the radial direction of the sphere itself. In 
the former case the gravitational potential comes from the Euclidean space in 
which the sphere is embedded, whereas in the latter case the intrinsic geometry 
of the surface affects the gravitational potential. 

A brief review of the literature relevant to eq. (|12|l — in particular, its relation 
to the nonlinear Hodge theory introduced in — is given in Sec. 1 of |10j . 
To those remarks we add that the extension of geometric variational problems 
for the Dirichlet energy to more general classes of energies was already outlined 
in 0], with a suggested application of the harmonic map energy to the theory 
of elasticity and that, whereas nonlinear Hodge theory generally considers the 
elliptic case of the variational equations, eq. I|12[l will be allowed to change from 
elliptic to hyperbolic type as Q cr it is exceeded. 

3 Conditions for trivial flow 

We will say that a flow is trivial if its flow potential u is a constant function. 
Equivalently, the velocity field of a trivial flow associates the zero vector to every 
point of the flow domain. 
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du 



We expect that a relative minimum for a smooth function of a single variable 
will occur at a point for which the second derivative is non-negative. The 
analogue of the second derivative for variational integrals is the second variation. 
Let u St t ■ M — > N, — e < s, T < e, be a compactly supported two-parameter 
variation such that ito.o = u. Define 

V = 

and 

W = 

The second variation is the quantity 



S.t I 



m 

du s j 
~ds~ 



\s,t=0 



\s.t=0- 



A map is said to be stable if /(V, V) is non-negative for any compactly supported 
vector field V along u. 

This definition implies the triviality of a certain class of flows: 

Theorem 1. Let the domain of a shallow, steady, irrotational hydrodynamic 
flow satisfying the hydrostatic law be represented by a compact Riemannian man- 
ifold M and let the flow potential u take M into the m-sphere § m for m > 2. Let 
the flow speed Q be given by 0) for < Q < 2. Then any stable flow potential 
takes every point of M to a single point on § m . 

Proof. We impose a geometric constraint on the variational problem (0)1,0, 
declaring that the image of the flow potential u must lie on a smooth, compact, 
Riemannian manifold N, where N is a submanifold of R fe for some sufficiently 
large number k. This transforms 10 into JHJ and the Euler-Lagrange equation 
(0J into the equation 



1 - | 1 du 



= [ 1- | )A(du. du). CI 4) 



where 8 g is the formal adjoint of the exterior derivative; g is the Riemannian 
metric on M; Q satisfies ©; A is the second fundamental form of N, where 
N is expressed as a submanifold of some higher-dimensional Euclidean space 
as noted earlier. The system (|14|) is identical to the variational equations for 
E p (u; M, N) with p given by (fH5|) in the special case C = n = 7 = 2, where C 
is the constant of Bernoulli's formula. 

Computing the second variation of E p (u; M, S m ) , we obtain 



d 2 f Q d r Q 

I\V.V)= I \du\'{\du\ 2 — I p(s)ds + (2-m)— J P (s)ds\dM. 



M 



dQ 2 



(15) 
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Because p(s) = 1 — s/2, the right-hand side of i|T5|) is negative under the hy- 
potheses of the theorem unless u is trivial almost everywhere. But u is contin- 
uous because \du\ is bounded, so the conclusion holds everywhere on M. This 
completes the proof. 

Remarks, i) The conditions on the domain M will only correspond locally 
to the physical model of Sec. 1.2. 

ii) Computing (|1(J|> for p given by l|13fl with C = n = 7 = 2, we find that 
Theorem 1 holds for flow speeds that extend well into the range of shooting 
flow. Moreover, solving Q = c 2 = 1 — Q/2 for Q, we verify that the Froude 
number attains the value 1 at Q = 2/3 = 2/(7+l) = Q cr it- 

in) A map u : M — ► N is said to be F -harmonic if it is a critical point of 
the -F-energy functional 



where F : [0, 00) — * [0, 00) is a strictly increasing, twice-differentiable function 
of its argument, inharmonic maps were introduced in [Tj as a unification of p- 
harmonic and exponentially harmonic maps. A version of inequality l|15|) holds 
for any map which is inharmonic; c.f. Theorem 7.1 of PP, taking 



In fact, the extension of Theorem 1 to a larger class of target manifolds is easily 
obtained by adapting the ideas of However, Theorem 1 is false for the best 
known special cases of i^-harmonic maps: p-harmonic, exponentially harmonic, 
and a-harmonic maps; c.f. the remark following Corollary 7.2 of [J. Moreover, 
critical points of (JHJl , © will not be inharmonic whenever the variational density 
e(u) is a decreasing function of Q. That possibility is allowed by condition 1)1 OJI 
and by our initial definition of p. 

iv) Theorem 1 recalls the famous statement that the wind can never be blow- 
ing simultaneously in the same direction everywhere on Earth. But of course 
the two assertions are mathematically different and have completely different 
proofs. They each result from the combination of a harsh global hypothesis — 
maximal symmetry in the range of a map — with the global imposition of what 
would be a reasonable local hypothesis — in one case, continuity and in the 
other, stability. 

v) Adding to the mass density p(s) a generalized "surface tension" of the 
form 



where p is a positive constant, does not affect the result of Theorem 1. But this 
extension is dependent on the sign of p. 





t{s) =/*(! + *) 



1/2 
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vi) Because the proof of Theorem 1 relies on applying a compressible model 
to shallow hydrodynamic flow, it can be viewed as a corollary of a theorem 
about compressible flow: 

Theorem 2. Let the compact Riemannian manifold M be the domain of 
a steady, polytropic, irrotational, perfect flow and let the flow potential u take 
M into the m-sphere § m for m > 2. Let the flow speed Q be given by 0) for 
< Q < 2/(7—1), where 7 is the adiabatic constant of the fluid. Then any 
stable flow takes every point of M into a single point ofS m . 

Proof. Follow the proof of Theorem 1 based on eq. (@J) (c.f. (2.14) of |3j), 
defining p as in (|13f) for any 7 > 1 and any n > 2. 

The correspondence between gas dynamics and shallow hydrodynamic flow, 
illustrated in 0} and leading to the similarity of Theorems 1 and 2, was ap- 
parently first reported in |KS| for time-dependent compressible flow in 1 space 
dimension. Elliptic-hyperbolic systems similar — or dual — to @ with a side 
condition of vanishing curl arise in projective geometry and optics, as well as in 
hydrodynamics and gas dynamics; c.f. 

This discussion provides a context for reinterpreting an earlier result. The 
variational arguments in [7j were given in terms of a function w (\du\ p ) . That 
function can be interpreted as the variational density of the map u, taking p = 2 
and 



Thus the results of |xj can be applied to hydrodynamic and compressible flow, 
which gives another set of conditions for triviality of the flow potential. 
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